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Recent studies exposed many remarkable properties of layered cobaltates Na^CoC^. Sur- 
prisingly, many-body effects have been found to increase at sodium-rich compositions of 
Na^CoCh where one expects a simple, nearly free motion of the dilute S — 1/2 holes doped 
into a band insulator NaCoC^. Here we discuss the origin of enigmatic correlations that turn 
a doped NaCoC>2 into a strongly correlated electronic system. A minimal model including 
orbital degeneracy is proposed and its predictions are discussed. The model is based on a 
key property of cobalt oxides - the spin-state quasidegeneracy of CoC>6 octahedral complex 
- which has been known, e.g., in the context of an unusual physics of LaCoOa compound. 
Another important ingredient of the model is the 90° Co-O-Co bonding in Na^CoCh which 
allows nearest-neighbor ti g — e g hopping. This hopping introduces a dynamical mixture of 
electronic configurations t| g ,S = and t% g e gi S — 1 of neighboring cobalt ions. We show 
that scattering of charge carriers on spin-state fluctuations suppresses their coherent motion 
and leads to the spin-polaron physics at x ~ 1. At larger doping when coherent fermionic 
bands are formed, the model predicts singlet superconductivity of extended s-wave symme- 
try. The presence of low-lying spin states of Co 3+ is essential for the pairing mechanism. 
Implications of the model for magnetic orderings are also discussed. 

§1. Introduction 

The physics of transition metal oxides offers a wealth of interesting phenomena 
related to their strongly correlated nature. This is because the bandwidth is rela- 
tively small compared to the intraionic Coulomb repulsion between the 3d electrons. 
The understanding of many unique properties of oxides, such as high-T c supercon- 
ductivity and colossal magnetoresistivity is therefore based on the well-known Mott 
physics. 1 ) 

Among the various families of transition metal oxides, attention has recently fo- 
cused on layered cobaltates, in particular on the sodium cobaltate Na x Co02- The in- 
terest was initiated by the large thermoelectric power observed in this compound, 2 )~ 4 ) 
i.e. the capability of an efficient conversion of heat energy to electricity. The research 
on these systems was further boosted by the unexpected discovery of superconductiv- 
ity (SC) in water-intercalated Na^OoC^. 5 ) Soon after, many remarkable properties 
were found 6 ) such as spin-sensitive thermopower, 3 ) unusual charge and spin order- 
ings, 7 )~ n ) very narrow quasiparticle bands 12 -* -16 ) and especially a very unusual phase 
diagram. 7 ) While the strongly correlated nature of Na^OoO^ is no longer at doubt, 
the mechanisms by which the correlated electrons design such an exotic phase dia- 
gram are not fully understood even on a qualitative level. 

Layered cobaltates consist of C0O2 planes with a triangular lattice of Co ions 



2 



Jin Chaloupka and Giniyat Khaliullin 



(a) z (b) 




Fig. 1. (a) Structure of a hexagonal C0O2 plane with 90° cobalt-oxygen bonds pointing along x, 
y and z directions. The directions in the hexagonal lattice of cobalt ions are denoted by a, 6 
and c respectively, (b) Each cobalt ion is surrounded by six oxygen ions forming an octahedron, 
in the real structure this is slightly distorted, (c) Each cobalt ion is a member of three C0O2 
chains which determine planes perpendicular to the x, y and z axes respectively. In the figure, 
chain of b direction determining a plane perpendicular to the y axis is shown. 

(see Fig. 1), separated either by Na layers as in Na x Co0 2 7) or by BiO-BaO lay- 
ers of rock-salt structure in so-called "misfit" cobaltates [see 16), 17) and references 
therein]. By suppressing electron motion along the c-axis these layers impose quasi- 
two-dimensionality upon the C0O2 layers. In addition they provide the charge car- 
riers to the C0O2 planes in doped compounds. Depending on the composition, the 
valence state of Co ions may vary in a wide range from non-magnetic Co 3+ t\ g S = 
state (as in NaCoC^) towards the magnetic Co 4+ t\ g 5 = 1/2 configuration (in 
Na x Co02 at small x). 

The Co 3+ systems where the cobalt ions have full t\g shell are naturally referred 
to as band insulators, 18 )' 19 ) whereas in Co 4+ 5 = 1/2 rich systems strongly correlated 
Mott behavior is expected. Thus by reducing sodium content in Na^CoC^ one should 
be able to observe an evolution from the weakly correlated band-insulator regime 
to a doped Mott limit. However, the experiments show completely opposite trends. 
The Co 4+ rich compounds are just moderately correlated metals 7 ) ,2 °) while Co 3+ 
rich compositions show signatures of strong correlations such as magnetic order, 7 )' 9 ) 
strong magnetic field effects, 3 ) etc. Pronounced incoherent structures in the angular- 
resolved photoemission (ARPES) spectra 16 )' 21 ) revealing the complex structure of 
holes doped into NaCo02 are also enhanced near the band- insulator limit. 

Even more puzzling situation occurs with hydration of Na x Co02. Depending on 
the amount of water intercalating the structure, Na x Co02 forms a monolayer hydrate 
or (at a larger water content) a bilayer hydrate 22 )' 23 ) (see Fig. 2 for the details of the 
structure). As a big surprise, superconductivity was discovered in the bilayer hydrate 
of Na a: Co02 with x « 0.35. 5 ) The corresponding transition temperature is quite low 
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N^CoOj - anhydrate monolayer hydrate bilayer hydrate 

Fig. 2. Layered structure of (a) anhydrate Na^CoC^ and water-intercalated (b) monolayer hydrate, 
(c) bilayer hydrate. Water-intercalation leads to an expansion in the direction normal to C0O2 
planes and in the case (c) brings a large screening of sodium potential. 

- T c ~ 5K in the optimal case. However, the identification of the pairing mechanism 
is a problem of principal importance, because this may shed some light on the other 
mysteries of Na^CoC^ as well. Soon after the discovery, comparisons with the cuprate 
situation have been made. 5 )' 24 ) They emphasize the quasi-two-dimensional C0O2 
layers of S = 1/2 background (Co 4+ ions) doped by S = (Co 3+ ions) charge 
carriers, thereby representing a ti g analog of CuC"2 planes. The triangular lattice 
of Co ions might be convenient for a realization of resonating-valence-bond (RVB) 
state, 25 ) which was originally proposed for a S = 1/2 Heisenberg antiferromagnet on 
a triangular lattice. 26 ) Despite this, it was quickly realized, that the phase diagram 
of Na^CoC^ 7 ) is radically different from that of cuprates. As already mentioned, 
the enhanced strongly correlated nature is rather observed at high relative content 
of nonmagnetic Co 3+ ions. Also the location of the superconducting dome initially 
expected at cobalt valency around 3.65 (as in Nan^CoC^) was later found at actual 
valency ~ 3.4, closer to the Co 3+ situation due to water-induced valency shift. 27 )~ 29 ) 

Many-body effects observed near a band insulator regime of cobaltates are puz- 
zling and indicate that electronic correlations originate from a specific mechanism 
different from that of the cuprates. We have recently proposed a model 30 )' 31 ) for 
strong correlations that operate over the entire phase diagram of Na^OoC^ includ- 
ing a weakly doped band insulator limit. In this paper, we discuss this model and its 
predictions in more detail, and further extend it to the case of multiorbital situation. 

Considering the limit of diluted 5 = 1/2 Co 4+ holes doped into non- magnetic 
band insulators NaCo02 and misfit cobaltates, the model predicts a polaronic behav- 
ior of the holes characterized by strongly reduced quasiparticle peaks and dispersive 
incoherent structures in their spectral functions. In the Fermi-liquid regime at large 
concentration of the doped holes, a specific pairing mechanism driven by correlated 
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Fig. 3. (a) 3d levels of cobalt ions split into fe 9 and e 9 levels in the octahedral field of oxygen ions. 
Trigonal distortion leads to a further splitting among the t^g orbitals. The order of the resulting 
ai s and e' g levels shown corresponds to that suggested by ARPES experiments in Na^CoC^. (b) 
Low-spin, intermediate-spin, and high-spin states of Co 3+ ions. The competition between the 
t'2 g — e g crystal field splitting favoring t 2g configuration and the Hund coupling favoring higher 
spin values makes them quasidegenerate. 



hopping of doped holes emerges in the model. The resulting SC is optimized near 
the cobalt valency 3.4 in agreement with experimental observations. The model is 
based on a unique aspect of Co 3+ ions, i.e., their spin-state quasidegeneracy, and on 
special lattice geometry of C0O2 layers realized in the layered cobaltates. 

Let us first consider the 3d electrons of cobalt ions from a single electron point 
of view as shown in Fig. 3(a). The individual cobalt ions are in octahedral crystal 
field [Fig. 1(b)] which leads to the usual t2 g — e g splitting WDq like in a perovskite 
structure. Due to the layered structure of cobaltates, which imposes a quasi-two- 
dimensionality of the electronic states, the symmetry of the originally three- fold 
degenerate t2 g levels is thus lowered and they split into so-called a\ g and e' g states. 
Crystal field prefers the t\ g configuration on the Co 4+ ions. The holes doped into 
NaCo02 or misfits are then expected to occupy the highest t2 g level. Initially, the- 
oretical considerations have suggested that this level should be of e' g symmetry (see 
Ref. 32) for details). However, experiments support an opposite scenario where the 
doped holes prefer to occupy the a\ g orbital states. Indeed, ARPES data shows a 
very simple Fermi surface derived from a single band of a± g symmetry. 16 )> 21 ) 

In the case of Co 3+ , the t% g S = configuration would be selected by the 
cubic crystal-field splitting IQDq alone. However, many-body effects represented by 
the Hund coupling favoring high-spin configurations compete with the crystal field 
splitting, making the S = 0, S = 1 and S = 2 states of Co 3+ quasidegenerate 
[Fig. 3(b)]. Spin-state quasidegeneracy of cobalt ions is well known, LaCoOs being 
a textbook example. 33 ) According to Ref. 34), magnetic states are in the range of 
~ 200 - 400 meV (S = 1) and > 50 meV (S = 2) above the ^5 = ground 
state (without lattice relaxations). This leads to a distinct property of Co 3+ rich 
oxides called "spin-state-transition" responsible for many anomalies such as spin- 
state changes in LaCoOs, 34 )' 35 ) and a spin-blockade effect in HoBaC^Os.s, 36 ) to 
mention a few manifestations of the Janus- like behavior of Co 3+ . 
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In the groundstate the cobalt ions are supposed to be in S = 1/2 Co 4+ or 
S = Co 3+ states. Thanks to the spin-state quasidegeneracy, a new, third degree 
of freedom - namely the low-lying S = 1 t^e^ state of Co 3+ - could be in principle 
employed by intersite t2 g — e g hopping which mixes the t\ g and t^e^ configurations 
of neighboring ions.*-* Whether such a matrix element is finite or not depends on 
lattice geometry. With the 180° Co-O-Co bond geometry of perovskite cobaltates, 
the t2 9 and e g sectors are separated with respect to the nearest neighbor hopping [see 
Fig. 4(a) for an explanation], which is therefore incapable to produce the S = 1 Co 3+ 
states. New situation encountered in layered cobaltates is that Co06 octahedra are 
edge-shared. In this geometry, the hopping occurs along the 90° Co-O-Co bonds, 
where the largest matrix element is that between the orbitals of t^g and e g symmetry 
as shown in Fig. 4(b). Doped holes can now easily generate S = 1 t^e^ states of 
Co 3+ which become strongly coupled to the groundstate by virtue of intersite t2 g — e g 
electron transfer. In other words, the magnetic configuration of Co 3+ is activated 
once the mobile Co 4+ holes are added in Na x Co02. A dynamical generation of t^e^ 
S = 1 states by a hole motion converts it into a many-body correlated object - the 
spin-polaron. At larger density of Co 4+ when spin-polarons overlap, the virtual 5 = 1 
states act as mediators of an effective spin-sensitive interaction. We eliminate these 

The S — 2 Co 3+ configuration with two e g electrons is not accessible by the ti g — e g hopping. 
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Fig. 4. (a) Hopping between d-orbitals of cobalt ions via oxygen ions in the case of 180° Co-O-Co 
bonds. The geometry allows for the nearest- neighbor ti g — t2 g and e g — e g hopping. The tz g 
and e g sectors are independent, (b) The same for the 90° Co-O-Co bonds as in Na^Co02- 
Nearest-neighbor t2 g — t2 g and t2 g — e g hopping is possible, the latter one being approximately 
two times stronger, (c) a, b and c directions in the hexagonal lattice of cobalt ions and orbitals 
connected by the t2 g — e g hopping along the three directions. 
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virtual states perturbatively, and find an effective model in a form of spin-selective 
pair hopping of electrons. The correlated hopping energy is optimized when holes 
are paired and condense into a SC state. 

A resume is that low- lying magnetic states of Co 3+ , accessible for electrons via 
the intersite hopping, provide an extra dimension in physics of Na^CoO^. In Sec. 2, 
we design a model incorporating this idea. Based on this model, we demonstrate in 
Sec. 3 that holes doped into the band insulator NaCo02 behave in fact as magnetic 
polarons dressed by the spin-state fluctuations of Co 3+ ions that are excited by 
the hole motion. Sec. 4 derives the effective interaction between holes, mediated 
by virtual spin-state excitations of Co 3+ ions, in a Fermi-liquid regime at larger 
hole densities. In Sec. 5, we discuss the relevance of these interactions to the spin 
ordering, and, by exploring their effect on the spin susceptibility, we find signatures 
of 2/cp-instabilities. Finally, we focus in Sec. 6 and 7 on the superconductivity within 
our model. First, in Sec. 6 we discuss symmetry and doping dependencies of pairing 
instabilities in the model based on the a\ g band only. To address a possible role of 
e' g pockets in the superconductivity, in Sec. 7 we extend the model by employing the 
full orbital structure of the relevant states. Sec. 8 concludes the paper. 

§2. Model Hamiltonian 

The t2 9 orbitals in Na^OoO^ split into single a\ g = (d xy + d yz + d zx )/y/3 state 
and two e' g states: e' gl = (d yz — d zx )j\[2 and e' g2 = (2d xy — d yz — d zx )/\^6. The 

photoemission experiments 12 ^ 14 ^' show that a single band, derived mostly from 
the aig orbitals, is active near the Fermi level [see Ref. 37) for possible orbital- 
selection mechanism]. Therefore, we consider first a model based on the a\ g = f 
hole states (the version containing also e' g orbitals will be introduced in Sec. 7). 
Valence fluctuations d^d^ — ► djdf (see the left part of Fig. 5) within the low-spin t 2ff 
manifold are described by the tight-binding Hamiltonian 

Ht = -tY,flfi*> (2-1) 

ija 

where t = 2i /3 and to = ^WA^ ^ s over l a P between i 2 g orbitals 38 ) (hereafter, a 
hole representation is used). The fermionic operators / are subject to a conventional 
Gutzwiller constraint (at most one hole at a given site). As a result the bandwidth 
following from (2-1) is reduced by the Gutzwiller factor g t = In&jiX + n^) (« at 
the relevant dopings). Here is the relative fraction of Co 3+ ions. 

Our crucial observation following from Fig. 4(b) is that the t2 g -e g hopping t = 
tfjtn/Apd, which uses a stronger a-bonding path with t a /t n ~ 2, leads to more 
effective valence fluctuations. These will constitute the second part, H^, of our 
model Hamiltonian: H t _i = H t + H^. The t process (Fig. 5) generates S = 1 state of 
Co 3+ composed of a t2 g hole and an e g electron. Low-spin S = t2g e g state is much 
higher in energy and can be ignored (see Sec. 7 for details on the t\ g e l g multiplet 
structure of Co ions). We represent the resulting state by the operator T specified 
by the spin projection of the 5 = 1 state and the e g orbital 7 created by t hopping, 
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Fig. 5. Processes contained in the t — i model. Usual i-hopping gives rise to the band structure, 
t-hopping allows to employ the low-lying triplet state of Co 3+ . High-spin (S = 2) state is not 
accessible by hopping. 



«-e., Tj 1>7 = e^/j, Tj lj7 = e^/j and T | 7 = (ej^/j + The e g orbital is 

selected by the hopping geometry depicted in Fig. 4(b). The nearest neighbor Co 
ions and two O ions binding them determine a plane [see Fig. 1(c)] which is labeled 
as a, b or c according to the Co-Co bond direction. With respect to this plane, the 
i-hopping couples the in-plane t2 g orbital to the out-of-plane e g orbital. This rule is 
shown in Fig. 4(c). For the t hopping along a, b, or c bond, 7 = 3x 2 — r 2 , 3y 2 — r 2 or 
3z 2 — r 2 respectively. The part of our minimal model Hamiltonian for Na x Co02 
reads then as 30 ) 

v ij 

(2-2) 

The factor of 1/V~3 comes from the projection onto the a\ g band within t2 g sector. 
Hi moves an electron from Co 3+ to Co^ + - producing a t2 g hole on site j - and 
replaces the t2 g hole on site i by a complex excitation T . Making use of the t2 g -e g 
hopping (the largest one for 90° Co-O-Co bonds), an electron "picks- up" the spin 
correlations present in virtual S = 1 states. 

The index 7 of the T excitation is determined by the orientation of the (ij) bond 
according to the rules in Fig. 4(c). The overlap between e g orbitals specified by 7 
and 7' is (7(7') = cos(^ 7 — Consequently, the excitations Try inherit the same 
overlap: (1^T\) oc (7|7'). This brings peculiar geometrical factors into the theory, 
as will be observed in the following sections. 

The T-excitation energy Ex is determined by all the many-body interactions 
within the C0O6 complex (Hund's coupling, p—d covalency, crystal field, . . . ) 34 ) This 
is a free parameter of the model. Experimentally, S = 1 states of CoOg complex 
in perovskite compound LaCoOs are found at energies E T ~ 0.2 - 0.4 eV 34 ) as 
already mentioned in the introduction. Based on this observation, we will use in this 
paper the representative value Et — 0.2 — 0.3 eV for layered cobaltates. In units 
of the ai g hopping integral t ~ 0.1 eV (which follows from the band structure fit 
to — 0.15 eV 37 - 1 ), this translates into Ex/t = 2 — 3 adopted below in our numerical 
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data. In principle, we expect some material dependence of Et as it is decided by the 
balance of several competing interactions. It is therefore highly desirable to quantify 
a multiplet structure of the CoOg-complex in Na^CoO^ as done in LaCo03. 34 -* For 
the ratio of the hopping amplitudes t and to, we set t/to = 2 as ta/t^ ~ 2. 

§3. Spin-state polaron behavior of quasiparticles 

Based on the model introduced in the previous section, we develop now a theory 
for the photoemission spectra in cobaltates at large sodium content. The concentra- 
tion of the Co 4+ holes doped into originally non-magnetic band insulator NaCoO-2 is 
small, so that we consider here the individual motion of the holes only. It is evident 
from the Hamiltonian (2-2), that by creating and destroying T excitations as they 
propagate, the holes are strongly renormalized and we deal with a spin-polaron prob- 
lem. This resembles the problem of doped Mott insulators like cuprates, however, 
the nature of spin excitations is different here because of the nonmagnetic ground 
state. Instead of magnon-like propagating modes as in cuprates, fluctuations of the 
very spin value of Co 3+ ions are the cause of the spin-polaron physics in cobaltates. 
In contrast to Refs. 39), 40), where a static hole surrounded by S = 1 Co 3+ ions was 
studied, in the present model the triplet S = 1 excitations are virtual and generated 
dynamically by the very motion of the hole via the t process. These two pictures 
can merge if a hole is strongly trapped (e.g., by Na-potential 41 )). 

For the calculation of the fermionic self-energies, we employ the self consistent 
Born approximation, which has extensively been used in the context of spin-polarons 
in cuprates. 42 ) First, we focus on spin excitation spectrum. Since a direct e g -e g 
hopping in case of 90°-bonds is not allowed by symmetry, the bare T spin excitation 
is a purely local mode, at the energy Et- The coupling to the holes in (2-2) shifts 
and broadens this level. Accounting for this effect perturbatively [see Fig. 6(a) 
for a diagrammatic representation], we obtain the T Green's function P _1 (io>) = 
iu — Et — Et(iuj) with 

2i 2 

E T (jw) = JpYl r k Qo{k,ie)Qo{k',ie + \uo) . (3-1) 

fcfc' ,ie 

Here, Go is the bare electron propagator Qo{k, ie) = (ie— ^fc)™ 1 with the a\ g dispersion 

on a triangular lattice £fc = —2t(c a + q, + c c ) + fi, where Cry COS kry cHlCl fvry are the 

projections of k on the directions a, b, c. The underlying E g symmetry of T operators 
involved in t hopping results in the factor = c 2 , + c 2 + c\ — c a Cb — CbC c — c c c a . We 
neglected a weak momentum dependence of Ut for the sake of simplicity. This is 
justified as long as St is small compared to the spin gap Et- 

Further, we approximate by its Brillouin-zone average 3/2, obtaining the 
simple expressions for Et in terms of bare fermionic density of states No(x) = 
E fc <^-6c): 

lmE T (E) = -Trt 2 / dx N {x)N (x + E) , (3-2) 

J-E 



Unusual electron correlations in Na x Co02 due to 



9 




Fig. 6. (a) Selfenergy of the T excitation in the lowest order approximation. Solid lines represent 
bare hole propagator, (b) Selfenergy of the holes in a self-consistent Born approximation. Double 
straight line is the full hole propagator. Double wiggly line represents the propagator of the T 
excitation renormalized by the selfenergy in (a). (After Ref. 31).) 



J-oo Jx 2 y - E 

These equations determine the renormalized spectral function of the spin-triplet ex- 
citation: pt{E) = — ■K~ 1 lmT>{\uj — > E+\8). Next, we use pt{E) below for calculation 
of the fermionic selfenergy. 

Holes are renormalized by creating and destroying S = 1 Co 3+ excitations while 
moving in the predominantly Co 3+ background. The selfenergy diagram in a self- 
consistent Born approximation accounting for this effect [Fig. 6(b)] reads as: 

£fc(ie) = --r E i r k>V(-™) + r h V(h>)] G(k',ie + iu) , (3-4) 

k ,io) 

where Q~ 1 (k,ie) = ie — — Zfc(ie). We can decompose the fc-dependence of the 
selfenergy into two simple terms and write Sk{oj) = 2t 2 [<P(uj) + i~fcH(u;)], where 



oo roo 



<P(u,)=l &Ep T {E)l dx * {X) , (3-5) 

Jo Jo to- E-x + id 

S(w) = HdEpriE) f dx ^ {X \.. . (3-6) 
Jo J-oo ^ + £ - X + 1(5 

The full local density of states N{E) = ^ fc A(k, E) and its E g symmetry part 
N(E) = Y^kEkA(k, E) are functions of the selfenergy itself, via the spectral func- 
tions A(k, E) = — -K~ l ImG(ie —>■ E + iS). The above equations are thus to be solved 
self-consistently. 

The reliability of the approximations made was tested by the comparison with an 
exact diagonalization of H t _^. Considering a single hole injected on a 7-site hexagonal 
cluster, it was shown in Ref. 31), that the above equations give results consistent 
with the exact diagonalization, even at rather small spin gap values Et ~ t. 

To illustrate the gross features of the hole renormalization, in Fig. 7 we show a 
complete map of the spectral function along M — T — K path in the Brillouin zone 
and spectral profiles at selected points. We have used a representative value Et = 2t 
which is renormalized by holes to Et ~ 1.4i. Renormalization of the holes leads to 
spectral functions with a reduced quasiparticle (qp-) peak whose spectral weight is 
transferred to a pronounced hump structure. [A peculiar momentum dependence of 
the matrix elements Ifc (note that Ifc = o = 0) reduces the effect at k = r point]. 
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Fig. 7. (a) Intensity map of the spectral function A(k, E) of the Co 4+ holes along the M — T — K 
path in the hexagonal Brillouin zone (inset) calculated at 30% doping and Et = It. Near 
the Fermi level, the quasiparticle (qp-) peak is well developed. As the hole energy reaches 
the renormalized T excitation energy Et, the quasiparticle peak broadens and its weight is 
transferred to a broad, incoherent background. This results in a "peak-dip-hump" profile of the 
spectral function as seen in the right panel. The top of the smoothed hump structure is indicated 
by triangles. The dashed line shows the bare dispersion, (b) Spectral function profiles at the 
r, K and M points and the total density of states N(E). The incoherent structure dominates 
N(E), a small coherent part resembling the density of states of the bare band (with a reduced 
bandwidth) is indicated by the gray area. Several maxima on the hump structure reflect the 
presence of multiple T excitations as sketched in the cartoon figures. (After Ref. 31).) 

Compared to the bare dispersion, the bandwidth of the renormalized holes is reduced 
by a factor ~ 2. The main observation here is that as the hole energy reaches Et, 
the dynamical generation of S = 1 excitations becomes very intense and a broad 
incoherent response develops, leading to the pronounced "peak-dip-hump" structure 
of A(k,E). Several maxima on the hump reflect the presence of multiple triplet 
excitations created by the hole propagation. All these are the typical signatures 
of polaron physics. Multiplet structure of the hump will in reality be smeared by 
phonons which are naturally coupled to the t transition involving also the orbital 
sector. Although experiments 43 ) indicate that electron-phonon coupling is moderate 
in cobaltates, it may enhance the spin-polaron effects as in cuprates. 44 ** Following 
the maximum of the smoothed hump structure, we observe its strong dispersion 
(stemming also from incoherent t hopping). 

Fig. 7 suggests a possible determination of Et from the quasiparticle damping. 
To address this problem, in Fig. 8 we show the energy width of the qp-peak following 
its dispersion curve. The sharp onset of the damping at the binding energy ps —lAt 
is clearly related to the maximum of the spin-excitation spectral function pt (E) . In 
addition, Fig. 8 shows the weight of the qp-peak which is fc-dependent (mainly due 
to the matrix element /&). 

Comparison of Figs. 7, 8 with the data of Refs. 16), 21) reveals a remarkable 
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Fig. 8. Top panel: Full energy width at half maximum of the quasiparticle peak along the M — F 
dispersion curve (see bottom-right panel) plotted as a function of the binding energy. The 
part from to M is accessible by ARPES experiments. Strong quasiparticle damping below 
~ —lAt is due to a scattering on S = 1 excitations. The spectral weight of the quasiparticle peak 
obtained by a direct integration is indicated by squares. When damping is small, it coincides 
with a conventional quasiparticle residue (1 — dE'/duj)^ 1 . Lower panel: The T-exciton spectral 
function (left), and renormalized hole dispersion (right). (After Ref. 31).) 

correspondence between theory and experiment. In particular, both the qp-peak 
and the hump dispersions (see Figs. 2, 3 of Ref. 16)) are well reproduced by theory, 
considering t ~ O.leV suggested by the band structure fit. 37 ) The onset energy Et ~ 
lAt for the qp-damping (Fig. 8) is then ~ 0.14 eV, in agreement with experiment 
(see Fig. 2(c) of Ref. 16)). 

Physically, at large x ~ 1 limit, dilute polarons are readily trapped by a random 
potential of Na- vacancies, 41 )' 45 ) thus qp-peaks should be suppressed at low hole 
doping. When the binding is strong, physics is local and a polaron takes a form of 
hexagon-shaped S = 1/2 object where a hole is oscillating to optimize both t and t 
channels. 39 -* Our model provides a microscopic basis for spin-polarons introduced on 
experimental grounds 8 )' n ) and discussed in detail in Refs. 39), 40). When the density 
of polarons is increased (as x decreases), they start to overlap forming narrow bands. 
Eventually, the polaron picture breaks down and a correlated Fermi-liquid emerges 
when x is further reduced. This is the subject of the following section, where we 
develop a perturbative theory accounting for the interactions between the holes. 

§4. Effective interaction between t2 g holes: single ai g band case 

In the Fermi liquid regime, the Eliashb erg- formalism, where the phonon shake- 
up processes (triggered by T-exciton) can also be incorporated, would be the best 
strategy. However, there are delicate constraints to handle: a lattice site cannot 
be occupied by two holes or by a hole and T-excitation simultaneously. For the 
sake of simplicity, we derive an effective fermionic interaction in a second order 
perturbation theory in t by considering the local virtual process depicted in Fig. 9. 
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Fig. 9. (a) Virtual process leading to the Hp Hamiltonian in Eq. (4-1). t-2g electron of a Cof + 
ion moves to e g level of the nearest-neighbor Co* + ion (process 1) and then to the t2 3 level of 
the next Co^ + neighbor (process 2). This is depicted in (b) as a motion of the hole-pair. The 



intermediate state contains the low-lying S — 1 state of Co 



3+ 



( After Ref. 30). 



This way, all the constraints in the intermediate states are treated explicitly. Such 
a perturbative treatment is valid as long as a polaron binding energy Ef, (an energy 
gain due to the t process) is small compared to a bare bandwidth W (~ 9t in a 
triangular lattice). From a self-consistent Born approximation discussed above, we 
obtained single polaron binding energy Ef, ~ 2.6i (2.9t) for Et set to the value 3i 
(2t) used in this paper. Hence, Ef, ~ and we can integrate out virtual spin 

states perturbatively. 

In the virtual process of Fig. 9 we consider an initial configuration with a pair of 
Co 4+ holes and a5 = Co 3+ state at the neighboring site. First t-hopping (i.e. the 
action of H^) generates an excited S = 1 state surrounded by two Co 4+ holes. The 
energy of this virtual state with respect to the initial state equals Et- The virtual 
state gets deexcited by second i-hopping leaving the system again with a pair of 
holes and a neighboring Co 3+ state. Since the pair now has a changed position, the 
virtual process effectively corresponds to the motion of the hole pair. The kinetic 
energy gain associated with this pair motion may lead to the pairing instability as 
shown in Sec. 6. Summing up the contributions of all possible virtual states, we 
arrive at the following effective Hamiltonian: 

H r = E [ P ^S\ 3 S k] + P&T^.Tj . (4-1) 

(ijk) 

The Hamiltonian in Eq. (4-1) describes the motion of the spin-singlet Sij = (/if/jj. — 
fiif n )/V2 and spin-triplet Ty = {frfr-U'rf,. + fM/V2, fa fa} Co 4 +-Co 4 + 
pairs in a background of S = Co 3+ ions. Sites i ^ k are the nearest neighbors 
of site j. No-double-occupancy constraint on / is implied when using this effective 
Hamiltonian. The pair-hopping amplitudes read as 

i%k = \y^iski-<t>jk), Pi3k = \p?3h- ( 4 - 2 ) 

We introduced here a parameter V = i 2 /Et characterizing the interaction strength. 
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The angles 4> £ (27r/3, 47r/3, 0) are selected by the orientation of the bonds (ij) and 
(jk) as already explained in Sec. 2. 

Two important remarks should be made on the pair-hopping amplitudes Pijk- 
First, they contain a geometrical factor cos(<^>jj — cj>jk) which equals +1 for straight 
pair hopping and —1/2 if the hopping process changes the direction of the pair. This 
factor is explained in Fig. 10. It originates from the fact, that different e g orbitals 
participate in the t-hopping process on bonds of different directions. Second, the 
pair hopping amplitude of singlet pairs is three times bigger than that of the triplet 
pairs. This nontrivial result originates from a quantum interference between different 
realizations of the virtual process as indicated in Fig. 11. The S = 1 intermediate 
state is fully transparent for singlets which equally use all three S z = ±1,0 states. 
However, these states contribute with different signs in case of triplets, resulting in 
a "spin blockade" for their motion. 

Finally, the processes contained in the reduced model Hamiltonian Ht-p = 
Ht + Hp, Eqs. (2-1) and (4-1), are summarized in a concise way in Fig. 12. In the 



(a) 



(b) 



xy (3r -r 2 ^ 2 -i 2 } = -1/2 

r*r I xy (D » 




Fig. 10. The origin of the geometrical factor cos((j>ij — <j>jk) in Eq. (4-1) due to the overlap of the e g 
orbitals in the intermediate state, (a) Along the c-bond, the xy and 3z 2 — r 2 orbitals are coupled 
via f-hopping. If the virtual process of Fig. 9 proceeds along the c-bond exclusively, we get the 
geometrical factor cos(0) = 1. (b) The bond-directions for the two i-hoppings differ. In the case 
of a-bond, the yz and 3a; 2 — r 2 orbitals are coupled via f-hopping. By taking the overlap of the 
e g orbitals active on the two bonds, we get the geometrical factor cos(27r/3) = —1/2. 
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Fig. 11. (a) Singlet-pair motion via an intermediate state composed of S — 1 Co 3+ ion (T exciton) 
and triplet state of the two holes at sites i and k. The process equally uses all three S z = ±1,0 
states of Co 3+ . The relative amplitude resulting from spin algebra equals 3/2. (b) The same for 
triplet-pair motion. Here the amplitude is distributed into two channels - with the two holes 
in the intermediate state being in singlet (one S z state of Co 3+ employed) or triplet state (two 
S z states of Co 3+ employed). Destructive interference of the two channels makes triplet-pair 
hopping amplitude 3 times smaller than that of the singlet-pair. 
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Fig. 12. Cartoon representation of the processes contained in the t — P Hamiltonian (2-l) + (4-l): 
(a) t-hopping of ai 3 holes, (b) pair-hopping of a singlet pair of holes, and (c) pair-hopping of a 
triplet pair of holes. 



following sections we will use this model containing: 1. the usual t-hopping of holes 
[Fig. 12(a)] and 2. pair-hopping of hole pairs [Fig. 12(b,c)] as the lowest-order effect 
of t-hopping. 

§5. Spin susceptibility 

The pair-hopping interaction in Eq. (4T) may be alternatively represented in 
a form of density-density and spin-spin couplings, emphasizing its relevance also to 
the charge and spin orderings: 

H P = V'^2 cos(4>ij - <f)j k ) [rijUik - \sjS ik ] . (5-1) 

(ijk) 

Here nj = \ f} a fja and s j = ^E ff a fjafja are the usual on-site charge and 
spin-density operators, but and with i ^ k are the charge and spin densities 
residing on bonds, i.e., n ik = \ ^ fj a f k(T , s z ik = \ Y, a vfLfka- Thus ' the interaction 
acts between the local (on-site) and non-local (bond) operators which is a consequence 
of its three-site nature. 

In momentum space, Eq. (5-1) can be written as 

H P = 2V^[n- q n q -\s- q ~s q ] , (5-2) 

Q 

with the usual operators n q = \ Y,k,a fl+ q ,ofk,*> 4 = I Dfc,a ° 'fl+ q ,Jk,<r> and the 
momentum counterparts of the non-local operators h q = \ Ylka^k+q,kf\. +qa f kcT , 
S q = h Efc, CT aF k+q,kfl +qj J kj(7 - The formfactor F k , k = cos(k a + k' a ) + cos(A: b + k' b ) + 
cos(A; c + k' c ) — c a c' b — Cbc' a — Cbc' c — c c c' b — c c c' a — c a c' c , where c' a = cos k' a , originates from 
a peculiar bond-dependence of interactions in Eq. (5-1). It manifests again that the 
h q and s q operators correspond to the particle-hole excitations that modulate the 
charge and spin bonds, respectively. 

To illustrate this unusual, nonlocal nature of correlations we investigate the 
effect of the interaction on the spin susceptibility within the RPA approximation. 
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local S non-local S 

Fig. 13. Diagrammatic representation of RPA equations (5-4) for the spin susceptibilities involving 
the interaction (5-2) between local (s q ) and non-local (s q ) spin densities that reside on sites 
and bonds respectively. Bare and RPA enhanced susceptibilities are represented by empty and 
shaded bubbles, respectively. 

The bare spin susceptibility Xss is given by the formula 

(ow 1 >p tanh(^ fc /2T) - tanh(q +(? /2T) 

A diagrammatic representation of the RPA approximation in our case of local- 
nonlocal spin interaction is shown in Fig. 13. RPA enhanced susceptibility follows 
from the equations 

Xss = xfj + A xfhis + A xf~Xss , (54) 

Xss = XsJ Ax's] Xss A XsS Xss 1 

where we defined a coupling constant A = 2V/3 = 2t 2 /3Ep. Because of the non- 
local spin density involved in the interaction, the corresponding bare susceptibilities 
X s °s = xfs an d xfs are also employed. They are obtained by multiplying the 
summed terms in Eq. (5-3) by a factor of F k+qk (for X s °s) or ^k+q k (^ or Xss)- ^he 
resulting spin susceptibility reads 



y (0) 
_ Xss 
Xss — jj 



D=(l-A X fJ) 2 -A\^ x fs ] - (5-5) 



The no-double-occupancy constraint imposed on Ht-p handled on a Gutzwiller level 
leads to a rescaling £ = gt£ and V = gpV. Since the interaction term comes from 
the hopping, the Gutzwiller factors are equal g t = gp pa and may be partially 
canceled in terms such as Ax^ leaving us with an effectively rescaled to and T. 

In Fig. 14 we show a sample result at the relative fraction of Co 3+ ions n& = 
0.5. The bare spin susceptibility [Fig. 14(a)] is concentrated around the r point. 
When the interaction is switched on [Fig. 14(b)], the 2kp ring in the susceptibility 
is enhanced. This suggests the fermionic 2/cp-instabilities in a Fermi-liquid phase, 
consistent with a picture inferred from the experiment. 46 ** Interestingly, the RPA- 
spin susceptibility at = 0.5 is most enhanced near the M point, i.e. near the 
observed magnetic Bragg peak position, 10 ) rather than at K typical for the AF 
Heisenberg spin system. In order to study the spin ordering at nd = 0.5 more 
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Fig. 14. (a) Map of bare x" /*+> for n d = 0.5 (the average Co-valency 3.5) at T = 0.025£ and w — > 0. 
(b) Corresponding RPA-enhanced susceptibility calculated at P/Et = 3t. The interaction 
enhances the susceptibility at the 2&f ring which (at given density rid = 0.5) nearly matches 
the Brillouin zone boundary. (After Ref. 30).) 



quantitatively, one should take into account also the Na ordering 7 ) which breaks a 
hexagonal symmetry of the underlying Fermi-surface. 

The 2/cf antiferromagnetic correlations within our model are manifested in the 
temperature dependence of the nuclear spin relaxation rate T^ 1 ~ (^ ? Imx/ w )m0' 
According to the Korringa law, in the normal state 1/T\T should be temperature 
independent. This is the type of behavior experimentally observed in monolayer 
hydrate of Na^CoC^. However, in bilayer hydrates, which show a superconducting 
transition at T c , 1/T\T is enhanced at low temperatures. 47 )' 48 ) This enhancement is 
stronger as T c increases, eventually leading to a magnetic order. 

We characterize the effect of the interaction (5-1) on 1/T±T by the enhancement 
factor (F 2 ) 



Using RPA enhanced Xss of Eq. (5-5) it can be approximately expressed as a double 
Fermi-surface average 

(F 2 )^(^ 2 [l-ARe( X ( S-F k ^)] 2 ) , (5-7) 

where D(q,ui), x^(q,^) and Xss{q,w) are evaluated at q = k' — k, lo —> 0. 

The resulting enhancement factor at the relevant fraction of Co 3+ ~ 0.5 
is presented in Fig. 15. To address the T < T c regime, we have included the ap- 
proximate BCS gap A = 1.76T C tanh(1.76Y / T' c /T — 1) in our calculations. As the 
interaction strength A approaches its critical value for a magnetic ordering, the 
enhancement factor dramatically increases resembling the experimentally observed 
critical behavior. 



§6. Superconductivity due to the pair-hopping interaction 

Now, we consider the Hp Hamiltonian in the context of superconductivity. The 
kinetic energy coming from the pair-hopping interaction in Hp is optimized when 
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Fig. 15. Enhancement factor (F 2 ) at n d — 0.5, t' /t = —0.2 as a function of temperature. With the 
interaction switched off, (F 2 ) = 1. The interaction enhances 1/TiT near the superconducting 
transition temperature as the coupling constant increases. The critical coupling constant for a 
magnetic ordering in this case equals /l cr it « At. 



Co 4+ holes move pairwise. Eventually, this leads to their condensation into SC 
state. The difference from cuprates is that pairs are formed here not due to the 
superexchange (in cobaltates, J is small 49 )) but because of the gain in the kinetic 
energy associated with t hoppings. It is evident from Eqs. (4T), (4-2) that spin- 
singlet pairs gain much more kinetic energy than triplets. Alternatively, it can be 
said that the S = Co 3+ states move more coherently when the S = 1/2 background 
is in a singlet state. This result, explained in Fig. 11, is a consequence of quantum 
interference among the rich variety of spin states involved during the virtual process 
leading to Eq. (4-1). 

A mean-field BCS analysis*) of Eq. (4-1) shows that Hp supports either extended 
s-wave singlet SC with the gap function oc 7(fc) = y / 2/3(c a + q, + c c ), or doubly- 
degenerate spin-triplet p-wave pairing with ^y x ,y(k) = {(s a — sj,), (2s c — s a — S{,)/\/3}, 
where s a = sin k a . The ci-wave channel is repulsive, while the /-wave one is attractive 
but too weak in the physically reasonable doping range. We estimated the T c from 

1= E %^«anh^, (64) 

where V a is either V or V/3, and the corresponding formfactors are j(k) or j x ,y(k) 
for the singlet s-wave and triplet p-wave pairing, respectively. To account for 
the no-double-occupancy constraint, the fermionic dispersion as well as the pair- 
hopping amplitude are renormalized by the Gutzwiller factor 37 ) g t = 2n^/ (l + n^) as 
(£, V") = (<7 t £, gtV), where rid is the relative fraction of Co 3+ ions. In the momentum 
summation, we have introduced a cutoff equal to the excitation energy Et- 

More details will be given in Sec. 7 in the context of the full orbital version of our model. 
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Fig. 16. (a) T c in the extended s-wave and p-wave channels. The complete profile of the dominant, 
s-wave T c curve is shown in the left inset together with 7^ (in arbitrary units) on the Fermi 
surface. The dashed (f)Vc = 1-5) and dotted (/3Vc = 3) T c curves are calculated including 
nearest-neighbor Coulomb repulsion which reduces the pairing interaction V at large rid- Shaded 
regions indicate the observed competing orderings (including the spin-charge order at rid = 0.5). 
(b) Probability ratio p(rid) (see the text for its definition) renormalizing the pairing interaction 
at different values of nearest-neighbor Coulomb repulsion relative to the effective temperature 
1//3 oc bandwidth. The feature at = 1/3 for large Vc manifests a honeycomb- lattice formation 
where each Co 3+ (o) has the maximum possible number of neighboring Co 4+ -Co 4+ pairs (•-•)• 
Above n d — 2/3, Co 4+ holes can avoid each-other completely if Vc is sufficiently large. (After 
Ref. 30).) 

We solved Eq. (6-1) at V = 3t (corresponding to t = Et = St). In terms of the 
BCS-coupling constant, this translates into A = VN = VN ~ 1/3 considering the 
density of states N ~ 1/W ~ l/9t. Therefore, the present formulation in terms of an 
effective fermionic Hamiltonian (4-1) should give reasonable results. At larger values 
of V, we encounter a strong coupling regime where one should use the original model 
(2-2) instead and treat virtual spin states explicitly. This limit remains a challenging 
problem for future study. 

The resulting T c values from Eq. (6T) are presented in Fig. 16(a) as solid lines. 
As expected, the highest T c values are found in the singlet channel, increasing with 
Co 3+ density due to the formfactor effect, until SC disappears at = 1 limit. A 
weak triplet pairing is present thanks to its formfactor matching well the Fermi sur- 
face, but it is expected to be destroyed by {e.g. Na) disorder. (We should notice that 
these trends are based on the present mean-field decoupling which ignores collective 
spin fluctuations. One can speculate, for instance, that the triplet pairing may be 
supported by ferromagnetic fluctuations within the C0O2 planes observed 9 ) at large 
rid limit). 

As the SC pairing considered here is due to the pair-hopping, Coulomb repulsion 
between the holes will oppose it. This is not a big trouble at high density of Co 4+ 
spins (as they cannot avoid themselves) but becomes a severe issue in a spin-diluted 
regime at large n^, where Coulomb repulsion reduces the process described in Fig. 9 
hence the amplitude V. Instead, the formation of spatially separated spin-polarons 
(Fig. 7) is favored, and competing orderings take over, such as an in-plane ferromag- 
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netism induced by a residual interactions between spin-polarons. 39 )' 40 -* To include 
the effects related to the Coulomb repulsion in the Gutzwiller fashion, we use an addi- 
tional multiplicative factor reflecting the suppression of the probability p^^ of having 
the required Co 3+ — Co 4+ — Co fe + configuration. We have determined this probability 
using a classical Monte-Carlo simulation of hardcore particles with nearest-neighbor 
Coulomb repulsion Vc moving on a hexagonal lattice of 1024 sites. The simulations 
were performed at different "effective temperatures" 1/(3 imitating the kinetic energy 
(of the order of bandwidth) which competes with the Coulomb repulsion in the real 
system. Plotted in Fig. 16(b) is the probability ratio p{nd) = Pijk(Yc) / Pijk(Yc = 0) 
for several values of (3Vc- The corresponding T c curves calculated with V — > p{nd)V 
locate the SC-dome near the valence 3.4, in a remarkable correspondence with ex- 
periment. 27 )~ 29 ) (The reported Co-valences ~3.4, 27 ) ~3.3, 28 ) ~3.46 29 ) optimal for 
SC translate into nd = 0.6,0.7,0.54). 

Finally, our t — t model provides a clear hint on the role of water-intercalation 
needed for SC in Na^CoC^. Without water, a random Na-potential induces some 
amount of spin-polarons locally (the origin of "Curie- Weiss metal" behavior 7 )) which 
suppress the pairing among the remaining fermions the usual way. Once this poten- 
tial is screened-out by the water layers, an intrinsic ground state of C0O2 planes as 
in Fig. 16 is revealed. (This interpretation of the water effect is consistent with the 
absence of superconductivity in the monolayer hydrate of Na x Co02, where the water 
resides in the Na layers.) The remaining "enemy" of SC is the Coulomb repulsion 
which prevents the pairing of dilute Co + fermions and supports the formation of 
spin-polarons and magnetism instead. 39 ) More pronounced polaron physics (because 
of the presence of large S = 1 T-exciton and narrow bandwidth) explains why T c in 
cobaltates is low compared to cuprates. Another mechanism for the water effect is 
provided by the band-structure calculations 50 -* that indicate a substantial flattening 
of the a\ g band-top and a reduction of the band splitting when the water-layers are 
present. To study the former effect, we have included negative t' in our calculation. 
Due to the combined effect of better formfactor utilization in the s-wave channel and 
Fermi velocity reduction this enhances singlet pairing as demonstrated in Ref. 30). 

§7. Full orbital structure of the interaction and the role of e' g pockets 

Local density approximation (LDA) bandstructure calculations on Na x Co02 
predict the existence of e' g parts of the Fermi surface. 51 -*' 52 ) According to some 
theoretical explanations of the pairing mechanism in Na x Co02, these so called e' g 
pockets may even play a crucial role. 53 ) -55 ) In the ARPES experiments, however, 
no e' g pockets are observed in the Fermi surface, and the corresponding band of e' g 
symmetry is well below the Fermi level in Na x Co02 at the relevant dopings. The 
reasons for such a discrepancy between LDA predictions and ARPES experiments 
became a highly debated topic and are not yet fully understood. 32 ) Nevertheless, the 
e' y band is reported to be located close to the Fermi level in the water-intercalated 
superconducting Na x Co02, 13 ) which opens the question of its possible role in super- 
conductivity. 
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Fig. 17. (a) Virtual process leading to the effective interaction between the holes. Occupied tig 
hole orbitals are labeled by a, [5, and 7 respectively. The f-active orbitals a and 7 are selected 
by the orientation of the (ij) and (jk) bonds, (b) Multiplet structure of a Co 3+ ion in t2 9 e g 
configuration. Two lower triplet states 3 Ti and 3 T2 differ in the occupied e g orbital being in- 
plane or out-of-plane with respect to the orbital occupied by the t2 g hole. B and C denote the 
Racah parameters. 56 - 1 



The aim of this section is to investigate the implications of the proximity of 
e' g states to the Fermi level within our model. To this end, we first derive a gen- 
eralization of the effective interaction Hamiltonian (4-1) considering the full orbital 
structure of the initial, intermediate, and final states. Using such a generalized inter- 
action Hamiltonian, we can study the case of mixed a\ g /e' g topmost band to assess 
the effect of its e' g parts. As the model bandstructure, we use the one from Ref. 37), 
where the multiorbital Gutzwiller approximation is applied to a tight-binding fit to 
LDA bands. The main advantage of this approach in the context of our study is the 
possibility to easily shift the position of e' g band with respect to the Fermi level by 
tuning a single tight-binding parameter. 

Instead of the projected a\ g holes we work here with the holes in canonical t2 g 
orbitals d xy , d yz , and d zx . It is convenient to use the orbital angles introduced in 
Fig. 4(c) to specify the orbitals, e.g., (4> a )i annihilates a hole in d yz orbital at site i. 
The a\ g hole operator / used in the previous sections corresponds to the linear 
combination fj a = [(4> a )ja + {^b)ja + (</ ) c)jo-]/v / 3 in the new notation. 

To derive an effective Hamiltonian of the form (4-1), we consider the virtual 
process depicted in Fig. 17(a) which is an ai 9 -unprojected equivalent of the pro- 
cess in Fig. 9. The initial state \i) of the two holes at sites k and j is either 
singlet Skj{^if3) = {ik^fijl ~ IkiPji) / V% or one of the triplet states T^j(^j3) = 
{lk]fij]-,{lk]l3ji + lkiPji)/y/2, IkiPji}- The t 2g orbitals 7 and [3 of the two holes 
are now explicitly indicated. The 7 orbital is the t-active t2 g orbital on (jk) bond. 
The amplitude of the transition to the final state |/) [either Sij(af3) or one of 
Tij(a(3)] contributing to the effective Hamiltonian Tip is evaluated as (f\Hp\i) = 
— ^ virt (/|Wj|virt)(virt|7ij:|i)/£^virt , where the summation runs over all possible in- 
termediate states. These consist of a hole pair (a hole at site % and 7 hole at site 
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k) and the excited Co 3+ ion at site j having the t^eg configuration. The excitation 
energy depends on the spin and orbital combination of the ti g hole and e g electron 
of the Co 3+ according to the multiplet structure presented in Fig. 17(b). The t- 
hopping Hamiltonian TC^ used here provides the hopping between tig and e g levels 
on nearest-neighbor Co ions as in Fig. 4(b,c). The resulting Tip Hamiltonian with 
full orbital structure can be written as 
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where a and 7 are selected by the directions of the bonds (ij) and (jk) respectively. 
The cosine and sine factors come from the overlap of the t-active e g orbitals at 
the respective bonds and the e g orbitals participating in the virtual states. In the 
following, we neglect the contribution of the high-energy S = states of Co 3+ . The 
simpler version of the Hamiltonian in Eq. (7-1) as given by Eq. (4-1) is obtained by 
letting E3 Tl = Ez T2 = Et and projecting on the a\ g states. 

Next we study the pairing interaction contained in (7T). In Sec. 6, it was shown, 
that the singlet pairing channel leads to the SC dome similar to the experimentally 
observed one. We therefore concentrate on the singlet-hopping term in (7-1) and 
estimate T c along the same lines as in Sec. 6. The corresponding singlet operators 
are first projected onto the topmost band (the closest to /i) fu a = Uka(4>a) ka + 
Ukb(<frb)ka + Uk, c {<t>c)ka as the relevant one for a T c estimation. A BCS Hamiltonian 



#BCS = J2 V kk'fk-\f-klf-k'l^ 



(7-2) 



kk' 



is then extracted from the result. The formfactor V kk * has a complicated form 
containing the bandstructure coefficients u ka as well as the overlap factors in a 
triple sum over orbital angles 



V kk' = 5Z( 2 ~ $a~/)ukaUk!3U k < pU k ^ COS k a COS k^ 

a/37 



sin(a — (3) sin(7 — (3) + — - — cos(a — (3) cos(7 — (3) 



Es Tl 



Ez 



T 2 



(7-3) 



In principle, Eq. (7-3) in its present form can be used for the calculation of T c , 
but to get a general insight to the results, a symmetry analysis of this formfactor is 
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Table I. Character table for the C& v group and a basis of its irreducible representations derived 
from cosk a , smk a . These functions are orthonormal on the BZ, i.e., ~}2 k 7i(fe)7j(fe) = 5ij . 



representation 


E 


2C 6 


2C 3 


c 2 


3<7„ 


3a d 


function 




A 1 


1 


1 


1 


1 


1 


1 


7.W = 


rr 

J 3 (COS k a + COS fc 6 + COS fe c ) 


Si 


1 


-1 


1 


-1 


1 


-1 


7/(fe) = 


^/| (sin fc a + sin fcf, + sin k c ) 


Ej 


2 


1 


-1 


-2 








7 P i(fe) = 


sin fc a — sin kb 
















7p2(fe) = 


-j= (2 sin fc c — sin k a — sin kb) 


E2 


2 


-1 


-1 


2 








7di(fe) = 


cos fc a — cos kb 
















7d2(fe) = 


-j= (2 cos fc c — cos k a — COS kb) 



necessary. Guided by the group theory, we express Eq. (7-3) using symmetry-adapted 
functions. 

In Table I, we show the irreducible representations of the Cq v group being the 
symmetry group of the hexagonal lattice of Co ions. The factors cos k a and cos k'^ 
as well as the bandstructure coefficients u can be expressed in terms of functions 
belonging to A\ and E 2 representation of the C% v group {B\ and E\ is employed in 
the case of triplet pairing). The functions 7 S , 7^1, and 7^2 allow us to directly express 
cosfco- and cosfc^, the bandstructure coefficients u are converted to a symmetry- 
adapted form via the relations a\ g = (u a + ui, + u c )/\/~3, e' gl = (u a — itb)/V2, e' g2 = 
(2u c — u a — lib) I \/6 deduced from Table I. We then perform the summations over a, 
(3, and 7 and regroup*) the resulting terms to form products of symmetric functions 
of k and k! respectively. For example, for the A\ representation (corresponding 
to the extended s-wave symmetry) there exist four such functions ty s \ = j s ai g ai g , 
&s2 = 775 {ldia lg e' gl + ~f d2 a lg e g2 ), & s3 = ^ 7s {e' gl e' gl + e' g2 e' g2 ) , & s4 = ldie' gl e' g2 + 

2 ld2{e g i-e g2 ). 

Disentangled interaction formfactor splits into the sum of s-wave, di-wave and 
c?2-wave contributions, i.e., these channels are independent in our BCS Hamiltonian. 
The formfactor of the s-wave channel reads as 

y^-wave) = x^ T (k){M 1 + rM 2 )W{k') (7-4) 

where we have introduced the vector \P T (k) = [ 1 P s i(k), l P s2 (k), l P S 3(k),^ , S 4(k)], the 
coupling constant A = P/Es^, the ratio of the triplet excitation energies r = 



This procedure is quite involved, since Eq. (7-3) at this stage contains a huge number of terms. 
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Fig. 18. (a) Symmetry functions 7 s (fe), 7di(fe) and jd2(k) entering the singlet-channel BCS in- 
teraction (7-5) resulting from (7-1) in the simple oi 9 case, (b) Symmetry adapted functions 
entering the BCS interaction in the singlet, extended s-wave channel in the general case. The 
bandstructure from Ref. 37) at = 0.5 was used here. 



Ez Tl jEz T2 and the matrices 
/ -2 -y/2 +v / 2 

Mi = - 



-V2 
-V2 



-1 
+1 
-1 



+1 
-1 
+1 



-V2\ 
-1 
+1 
-1 / 



M> 



2 — T 





f " 2 


-3^2 


-V2 


+V2\ 


1 


-3y/2 


+1 


+7 


-7 


4 




+7 


+9 


-9 




\+V2 


-7 


-9 


+9 / 



Each of the two d-wave channels could be represented in a similar way using 6x6 
matrices and corresponding six symmetry adapted functions. It is instructive to 
see, how a simple a± g expression emerges from the general result. In such case 



dig 



1 and e' 



-92 



everywhere in the Brillouin zone implying W. 



si 



7s 



]p S 2 = = & S 4 = and a similar reduction in the ci-wave part. The simplified 
expression for the singlet-pairing formfactor could be written as 

v kk' = — 



-A [- 7s (fc) 7s (fc') + 27 d i(fc) 7d i(fc') + 2 ld2 (k) ld2 (k')] (7-5) 



which shows an attractive extended s-wave part and a repulsive ci-wave part. The 
symmetry functions 7s , 7 ^i, and 7c f 2 are presented in Fig. 18 along with the functions 
&si, &S3, and for the bandstructure of Ref. 37). 

Finally, we estimate T c in the singlet, extended s-wave channel using the dis- 
entangled formfactor (7-4). Plugging in our form of V kk ' into the gap equation at 
T = T c , we find that the transition temperature T c is determined by the condition 
of the largest eigenvalue of the 4x4 matrix 



-<7tA(Mi + rM 2 ) £ ^(fc)^(fc) tanh( o f /2rc) 



(7-6) 
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Fig. 19. (a) Transition temperature T c as a function of rid for the bandstructure of Ref. 37). The 
aig-e' s orbital splitting A was modified to shift the e' g bands closer to the Fermi level. E3 Tl = 
t/2, r = 1/3 was used in the calculation. The effect of the nearest-neighbor Coulomb repulsion 
is not included. The inset shows the gap function in the Brillouin zone, (b) The effect of the 
modified ai g -e' g orbital splitting on the bandstructure at n d — 0.55. Topmost band employed 
in T c calculation is drawn as a solid line, (c) Distance of the top of the e' g band to the Fermi 
level as a function of na- 

being equal to 1. Here is the dispersion of the renormalized topmost band. 

Shown in Fig. 19 is the doping-dependent transition temperature as found from 
(7-6) employing the topmost band of the bandstructure of Ref. 37). To study the 
effect of the proximity of e' g band to the Fermi level we have varied the a\ g -e' g 
orbital splitting determining the position of the e' g band. As observed in Fig. 19, T c 
in the relevant range increases as the e' g band approaches fj, very closely, to around 
10 - 30 meV. 

Another effect of the e' g band is a strong anisotropy of the superconducting 
gap. In the pure a\ g case, the gap function as well as the nearest neighbor tight- 
binding dispersion is proportional to 7 a (fe), so that the gap is totally isotropic at the 
Fermi surface. In the present case, the gap consists of the four symmetry functions 
presented in Fig. 18(b). Near the Brillouin zone center, where the topmost band is of 
dig character, the gap is again determined by ^y s (k). However, closer to the Brillouin 
zone boundary, the band switches to mainly e' g character and the functions \P S 2, & S 3 
and ^4 start to dominate the gap reversing its sign. They bring strong anisotropy to 
the gap function at the Fermi surface which attains 30 — 60% at = 0.55 depending 
on the aig-e' splitting. Such a strong anisotropy of the gap may lie behind the 
absence of the coherence peak in the nuclear spin relaxation rate data. 57 )' 58 ) 

§8. Conclusions 

We have discussed a strongly correlated model for Na x Co02 which is based on 
the spin-state quasidegeneracy of S = and 5=1 states of Co 3+ and the specific 
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geometry of the C0O2 planes in layered cobaltates. The basic idea of the model is 
that t2 9 — e g hopping of 3d electrons between nearest-neighbor cobalt ions, which is 
enabled thanks to the 90° geometry of Co-O-Co bonds, allows to employ the low- 
lying S = 1 t^gCg state of Co 3+ . This new degree of freedom is exploited by the 
Co 4+ holes as they move in the C0O2 plane. 

In the sodium-rich region, when the Co 4+ holes are dilute, the model naturally 
explains experimentally observed strong correlations and interprets them in terms 
of a spin-polaron formation. At higher concentration of the doped holes, when the 
Fermi-liquid regime is established, we have derived effective interactions between 
the holes and discussed their impact on spin fluctuations. The superconductivity 
mediated by the spin-state fluctuations of Co 3+ ions emerges in the model, at exper- 
imentally observed compositions. Finally, we discussed the symmetry of the effective 
interaction in the context of the possible role of e' g bands in the superconductivity. 

Given the simplicity and experimentally motivated design of the model, its suc- 
cess can hardly be accidental. Therefore, t — t Hamiltonian, Eqs. (2T) and (2-2), can 
be regarded as a basic minimal model for Na^OoO^. The idea of spin-state fluctu- 
ations may also have broader applications, e.g., in oxides of Rh and Ir ions with a 
similar spin-orbital structure and lattice geometry. 

Formally, systems like LaCoO"3 and NaCo02 fall into the category of band in- 
sulators because Co 3+ ions have an even number of electrons forming a completely 
filled, spinless configuration t\ g in their ground state. Nevertheless, they are strongly 
correlated materials due to a proximity of virtual magnetic configurations of Co 3+ 
which can easily be activated by doping, temperature, etc. We conclude that the 
presence of such low-lying magnetic states is responsible for the rich and nontrivial 
physics of cobaltates. 
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